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Abstract 

An inclusion of temperature and chemical potential dependent surface tension into 
the gas of quark-gluon bags model resolves a long standing problem of a unified descrip- 
tion of the first and second order phase transition with the cross-over. The suggested 
model has an exact analytical solution and allows one to rigorously study the vicin- 
ity of the critical endpoint of the deconfinement phase transition. It is found that 
at the curve of a zero surface tension coefficient there must exist the surface induced 
phase tranition of the 2"'^ or higher order. The present model predicts that the critical 
endpoint (CEP) of quantum chromodynamics is the tricritical endpoint. 
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1 Introduction 



The role of surface tension for the quark gluon plasma (QGP) was discussed long ago [1, 2], 
however, up to now the situation is somewhat unclear [3, 4]. In nuclear and cluster physics 
the importance of the surface tension for the properties of CEP is known from a number 
of exactly solvable cluster models with the 1** order phase transition (FT) which describe 
the critical point properties very well. These models arc built on the assumptions that the 
difference of the bulk part (or the volume dependent part) of free energy of two phases 
disappears at phase equilibrium and that, in addition, the difference of the surface part (or 
the surface tension) of free energy vanishes at the critical point. The most famous of them 
is the Fisher droplet model (FDM) [6, 7] which has been successfully used to analyze the 
condensation of a gaseous phase (droplets of all sizes) into a liquid. 

Another such a model is a simplified version of the statistical multifragmentation model 
(SMM) [8] which was solved analytically both for infinite [9, 10] and for finite [11, 12] volumes 
of the system. The analysis of critical indices of the SMM [10] shows that the value of Fisher 
exponent tsmm — 1.825 ±0.025 of this model is consistent with ISiS Collaboration data [13] 
and EOS Collaboration data [14]. Such an experimentally obtained range of the r index is 
of a principal importance because it gives a very strong evidence that the SMM, and, thus, 
the nuclear matter, has a tricritical cndpoint rather than a critical cndpoint [9, 10]. 

This success of the SMM initiated the studies of the surface partitions of large clusters 
within the Hills and Dales Model [15, 16] and led to a discovery of the origin of the tempera- 
ture independent surface entropy similar to the FDM. As a consequence, the surface tension 
coefficient of large clusters consisting of the discrete constituents should linearly depend on 
the temperature of the system [15] and must vanish at the critical cndpoint. However, the 
present formulation of the Hills and Dales Model [15, 16], which successfully estimates the 
upper and lower bounds of the surface deformations of the discrete physical clusters, does 
not look suitable for quark-gluon bags. Therefore, in this work I insert the surface tension 
into the gas of bags model (GBM) [17], assume a certain dependence of the surface tension 
coefficient on temperature and baryonic chemical potential. Then I analyze the quark gluon 
bags with surface tension (QGBST) model and concentrate on the impact of surface tension 
on the properties of the deconfinement phase diagram and the QCD critical endpoint. 

Here I show that at low values of the baryonic chemical potential the I''* order deconfine- 
ment PT degenerates into a cross-over, if the surface tension coefficient becomes negative 
for lower values of temperature than the transition temperature. Also I prove the existence 
of an additional PT of the 2"^^ or higher order along the curve where the surface tension 
coefficient vanishes. Thus, I am arguing that the QGBST model predicts the existence of 
the tricritical rather than critical endpoint. 
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2 The Role of Surface Tension at Zero Baryonic Den- 
sities 



The isobaric partition of the QGBST model obtained from the grand canonical one Z{V, T) 
is as follows 

oo 

Z{s,T) = !JVeM-sV) Z{V,T) = j^^^y (1) 

Here the function F{s, T) consists of two parts, the discrete mass- volume spectrum Fh{s, T), 
and the continuous part of the spectrum Fq{s,T) 

F{s,T) ^ mm,)e-^^' + u{T) J dv -p[(-'^(^)-;>-(^)-1 , (2) 

Vb 

where the function (j){T,mk) = ^ Jp dp e r = 9k^^ -^2 (^) is the particle 

density of bags of mass and eigen volume and degeneracy Qk- 

At the moment the particular choice of function Fq{s,T) in (2) is not important. The 
key point of my treatment is that it should have the form of Eq. (2) which has a singularity 
at s = Sq because for s < sq the integral over the bag volume v diverges at its upper limit. 
As will be shown below the isobaric partition (1) has two kind of singularities: the simple 
pole s = s*fj and the essential singularity s = sq The rightmost singularity defines the phase 
in which matter exists, whereas a PT occurs when two singularities coincide [17, 9, 3]. All 
singularities are defined by the equation 

s* = F{s*,T). (3) 

The T;-linear term in the exponential of the continuous spectrum (2) is nothing else, but 
a difference of the bulk free energy of a bag of volume v, i.e. —Tsv, which is under external 
pressure Ts, and the bulk free energy of the same bag filled with QGP, i.e. —Tsqv. The 
term —Tsv appears due to the hard core repulsion [3], whereas the QGP pressure, Tsq, 
appears in (2) as a generalization of the Hagedorn mass spectrum [18]. At phase equilibrium 
this difference of the bulk free energies vanishes and the properties of phase equilibrium are 
defined by the surface free energy. 

Note that the usage of the grand canonical description for the exponential mass or volume 
spectrum of Hagedorn type was strongly criticized recently [4, 5, 19, 20, 21] because of the 
thermostatic properties of this spectrum. However, the hard core repulsion compensates the 
growing part of the mass- volume spectrum and, hence, the criticism of Refs. [4, 5, 19, 20, 21] 
is irrelevant to the present model. 

The new element in (2) is the presence of surface free energy (Tqv'^ (k < 1) of the bag. The 
power K < 1 which describes the bag's effective surface is a constant which, in principle, can 
differ from the typical FDM and SMM value |, if the highly non-sperical bags are possible 
[3, 7, 15, 16]. The ratio of the temperature dependent surface tension coefficient to T (the 

reduced surface tension coefficient hereafter) which has the form a{T) — ^ 



T 

{k — 0, 1, 2, ...). Here Uo > can be a smooth function of the temperature, but for simp! 
I fix it to be a constant. 



2fc+l 

icity 
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In choosing such a simple surface energy parameterization I follow the original Fisher 
idea [6] which allows one to account for the surface energy by considering some mean bag of 
volume V and surface v'^. The consideration of the general mass- volume-surface bag spectrum 
is reserved for the future investigation. In contrast to the FDM and SMM, the power n < 1 
which describes the bag's effective surface is a constant which, as mentioned above, can 
differ from the typical FDM and SMM value |. This is so because near the deconfinement 
PT region the QGP has low density and, hence, like in the low density nuclear matter 
[22], the non-sperical bags (spaghetti-like or lasagna-like [22]) can be favorable (see a [3] 
and references therein). A similar idea of "polymerization" of gluonic quasiparticles was 
introduced recently [23]. 

The second essential difference with the FDM and SMM surface tension parameterizations 
is that the vanishing of (j{T) above the CEP temperature is not required. As will be shown 
later, this is the most important assumption which, in contrast to the GBM, allows one 
to naturally describe the cross-over from hadron gas to QGP. Note that negative value 
of the reduced surface tension coefficient cr(T) above the CEP does not mean anything 
wrong. As discussed above, the surface tension coefficient consists of energy and entropy 
parts which have opposite signs [6, 15, 16]. Therefore, cr{T) < does not mean that the 
surface energy changes the sign, but it rather means that the surface entropy, i.e. the 
logarithm of the degeneracy of bags of a fixed volume, simply exceeds their surface energy. 
In other words, the number of non-spherical bags of a fixed volume becomes so large that the 
Boltzmann exponent, which accounts for the energy "costs" of these bags, cannot suppress 
them anymore. 

Finally, the third essential difference with the FDM and SMM is that it is assumed 
that the surface tension in the QGBST model vanishes at some line in — T plane, i.e. 
Tcep = Tcepins)- However, in the subsequent sections I will consider Tcep = Const for 
simplicity, and in Sect. 4 I will discuss the necessary modifications of the model with 

Tcep — Tcepi^fJi b) ■ 

In principle, besides the bulk and surface parts of free energy, the continuous vohime 
spectrum Fq{s,T) could include the curvature part as well, which may be important for 
small hadronic bubbles [24, 25, 26] or for cosmological PT [27]. It is necessary to stress, 
however, that the critical properties of the present model are defined by the infinite bag, 
therefore the inclusion into the function Fq(s, T) of a curvature term of bag's free energy of 
any sign could affect the thermodynamic quantities of this model at s = sq(T) and cr{T) — 0, 
which is possible at (tri)critical endpoint only (see below). If, the curvature term was really 
important for the cluster models like the present one, then it should have been seen also 
at (tri)critical points of the FDM, SMM and many systems described by the FDM, but 
this is not the case [7, 28]. Indeed, recently the Complement method [28] was applied to 
the analysis of the largest, but still mesoscopic drop of a radius Rdr representing the liquid 
in equilibrium with its vapor. The method allows one to find out the concentrations of 
the vapor clusters in finite system in a whole range of temperatures and determine the free 
energy difference of two phases with high precision. The latter enables us not only to extract 
the critical temperature, surface tension coefficient and even the value of Fisher index r of 
the infinite system, but also such a delicate effects as the Gibbs-Thomson correction [29] to 
the free energy of a liquid drop. Note that the Gibbs-Thomson correction behaves as R^r^ 
but the Complement method [28] allows one to find it, whereas the curvature part of free 
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energy, which is proportional to Rdr, is not seen both for a drop and for smaller clusters. 
Such a result is directly related to the QGP bags because QCD is expected to be in the same 
universality class [30, 31] as the 3-dimcnsional Ising model whose clusters were analyzed in 
[28]. Therefore, admitting that for finite QGP bags the curvature effects may be essential, I 
leave them out because the critical behavior of the present model is defined by the properties 
of the infinite bag. On the other hand, similarly to the FDM, SMM and FDM-like systems 
discussed in Ref. [7], I assume that the curvature part of free energy of the infinite QGP 
bag is not important and leave for the future analysis the question why this is so. 

According to the general theorem [17] the analysis of PT existence of the GCP is now 
reduced to the analysis of the rightmost singularity of the isobaric partition (1). Depending 
on the sign of the reduced surface tension coefficient, there are three possibilities. 

(I) The first possibility corresponds to (t(T) > 0. Its treatment is very similar to the GBM 
parameterization of continuous spectrum with r > 2 [17]. In this case at low temperatures 
the QGP pressure Tsq{T) is negative and, therefore, the rightmost singularity is a simple pole 
of the isobaric partition s* = sh{T) = F{sh{T),T) > sq{T), which is mainly defined by a 
discrete part of the volume spectrum Fh{s, T). The last inequality provides the convergence 
of the volume integral in (2) (see Fig. 1). On the other hand at very high T the QGP 
pressure dominates and, hence, the rightmost singularity is the essential singularity of the 
isobaric partition s* = sq{T). The phase transition occurs, when the singularities coincide: 

s,(T,)^'^ = s,m^^-^, (4) 

which is nothing else, but the Gibbs criterion. The graphical solution of Eq. (3) for all these 
possibilities is shown in Fig. 1. Like in the GBM [17], the necessary condition for the PT 
existence is the finiteness of Fq{sq{T),T) at s = Sq{T). It can be shown that the sufficient 
conditions are the following inequalities: Fq{sq{T),T) > sq{T) for low temperatures and 
F(sq(T),T) < sq{T) for T — > oo. These conditions provide that at low T the rightmost 
singularity of the isobaric partition is a simple pole, whereas for hight T the essential sin- 
gularity sq{T) becomes its rightmost one (see Fig. 1 and a detailed analysis of the case 

f^B 0). 



The PT order can be found from the T-derivatives of sh{T). Thus, differentiating (3) 
one finds 

/ _ G + ulCr-i{A,-a) ■ s'q 
- 1 + ^/C._i(A,-a) ' 
where the functions G and /C,— a(A, —a) are defined as {A = sh — sq) 

G^F'„ + iFQ+ ^ ^-«( A, -a) , (6) 

JCr-aiA, -a) ^ 7dv e-p[-^;::(^)"-] . (7) 

Vo 

Now it is easy to see that the transition is of the 1** order, i.e. s'q{Tc) > s'h{T^, 
provided (t{T) > for any r. The 2"'^ or higher order phase transition takes place provided 
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Figure 1: Left panel. Graphical solution of Eq. (3) which corresponds to a PT. The solu- 
tion of Eq. (3) is shown by a filled hexagon. The function F{s,C,) is shown by a solid curve 
for a few values of the parameter ^. The function F{s,C,) diverges for s < sq{^) (shown 
by dashed lines), but is finite at s = sq{^) (shown by black circle). At low values of the 
parameter ^ = C,a, which can be either T or fiB, the simple pole sh is the rightmost singu- 
larity and it corresponds to hadronic phase. For = C,b ^ the rightmost singularity is an 
essential singularity s = sq{^b), which describes QGP. At intermediate value ^ = both 
singularities coincide sh{^c) = sgi^c) and this condition is a Gibbs criterion. 
Right panel. Graphical solution of Eq. (3) which corresponds to a cross-over. The no- 
tations are the same as in the left panel. Now the function F(s,^) diverges at s = sq(^) 
(shown by dashed lines). In this case the simple pole Sh is the rightmost singularity for any 
value of ^. 

s'q{Tc) = s'fj{Tc) at T = Tc. The latter condition is satisfied when JCr-i diverges to infinity 
at T ^ (Tc — 0), i.e. for T approaching from below. Like for the GBM choice (2), such a 
situation can exist for cr(Tc) = and | < r < 2. Studying the higher T-derivatives of sh{T) 
at Tc, one can show that for cr{T) = and for {n + l)/n < t < n/{n — 1) {n = 3,4,5, ...) 
there is a n*^ order phase transition 

snin) = SQin), s'b{T^) = s'q{T^\ ... st'\T^) = s^q-'\t^), s^n\T^) + s^-\t^), (8) 

with s^^\Tc) = oo for {n + l)/n < r < n/{n — 1) and with a finite value of s^^\Tc) for 
r = {n + l)/n. 

(II) The second possibility, cr(T) = 0, described in the preceding paragraph, does not give 
anything new compared to the GBM [17]. If the PT exists, then the graphical picture of 
singularities is basically similar to Fig. 1. The only difference is that, depending on the PT 
order, the derivatives of F(s, T) function with respect to s should diverge at s = sq(Tc). 

(III) A principally new possibility exists for T > Tcep, where cr(T) < 0. In this case there 
exists a cross-over, if for T < Tcep the rightmost singularity is sh{T), which corresponds to 
the leftmost curve in the left panel of Fig. 1. Under the latter, its existence can be shown 
as follows. Let us solve the equation for singularities (3) graphically (see the right panel of 
Fig. 1). For o"(T) < the function Fq{s,T) diverges at s = sq{T). On the other hand, the 
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partial derivatives — ^y' < and — < are always negative. Therefore, the function 
F{s,T) = Fh{s,T) + Fq{s,T) is a monotonically decreasing function of s, which vanishes 
at s — > cxo. Since the left hand side of Eq. (3) is a monotonically increasing function of 
s, then there can exist a single intersection s* of s and T) functions. Moreover, for 
finite sq{T) values this intersection can occur on the right hand side of the point s = sq{T), 
i.e. s* > sq{T) (see the right panel of Fig. 1). Thus, in this case the essential singularity 
s — sq{T) can become the rightmost one for infinite temperature only. In other words, the 
pressure of the pure QGP can be reached at infinite T, whereas for finite T the hadronic 
mass spectrum gives a non-zero contribution into all thermodynamic functions. Note that 
such a behavior is typical for the lattice QCD data at zero baryonic chemical potential [32]. 

In terms of the present model it is clear that a cross-over existence means a fast transition 
of energy or entropy density in a narrow T region from a dominance of the discrete mass- 
volume spectrum of light hadrons to a dominance of the continuous spectrum of heavy 
QGP bags. This is exactly the case for (j{T) < because in the right vicinity of the 
point s = sq{T) the function F{s,T) decreases very fast and then it gradually decreases 
as function of s-variable. Since, Fq{s,T) changes fast from F{s,T) ~ Fq{s,T) ~ sq{T) to 
F(s,T) ~ Fh(s,T) ~ sh{T), their s-derivatives should change fast as well. Now, recaUing 
that the change from F[s,T) ~ Fq{s,T) behavior to F[s,T) ~ Fh{s,T) in s-variable 
corresponds to the cooling of the system (see the right panel of Fig. 1), 1 conclude that that 
there exists a narrow region of temperatures, where the T derivative of system pressure, 
i.e. the entropy density, drops down from ^ ~ sq{T) + T^^^ to || - sh{T) + T^^jf^ 
very fast compared to other regions of T, if system cools. If, however, in the vicinity of 
T = Tcep — the rightmost singularity is sq(T), then for T > Tcep the situation is different 
and the cross-over does not exist. A detailed analysis of this situation is given in Sect. 4. 

Note also that all these nice properties would vanish, if the reduced surface tension 
coefficient is zero or positive above Tcep. This is one of the crucial points of the present 
model which puts forward certain doubts about the vanishing of the reduced surface tension 
coefficient in the FDM and SMM. These doubts are also supported by the first principle 
results obtained by the Hills and Dales Model [15, 16], because the surface entropy simply 
counts the degeneracy of a cluster of a fixed volume and it does not physically affect the 
surface energy of this cluster. 



3 Generalization to Non-Zero Baryonic Densities 

The possibihties (I)-(III) discussed in the preceding section remain unchanged for non-zero 
baryonic numbers. The latter should be included into consideration to make our model more 
realistic. To keep the presentation simple, I do not consider strangeness. The inclusion of 
the baryonic charge of the quark-gluon bags does not change the two types of singularities 
of the isobaric partition (1) and the corresponding equation for them (3), but it leads to the 
following modifications of the Fh and Fq functions: 

Fh {s, T, i^b) = gje^-^^'HT, m,) , (9) 
Fq {s, T, iib) = u{T, hb) J dv — ^ ^ ^ . (10) 
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Here the baryonic chemical potential is denoted as /^s, the baryonic charge of the j-th hadron 
in the discrete part of the spectrum is bj. The continuous part of the spectrum, Fq can be 
obtained from some spectrum p{m, v, b) in the spirit of Ref. [33], but this is not the aim of 

the present work. 

The QGP pressure pg = Tsq{T, /ib) can be also chosen in several ways. Here 1 use the 
bag model pressure 



PQ = fo^ 



95 I 10 

T + ^ 



f /£b\^ _|_ 5 !iB_\ 
\ T J ^ 9^ \ T j 



(11) 



but the more complicated model pressures, even with the PT of other kind like the transition 
between the color superconducting QGP and the usual QGP, can be, in principle, used. 
It can be shown [3] that the sufficient conditions for a PT existence are 



F{{sQ{T,fiB = 0) + 0),T,fiB = 0) > SQ{T,fiB = 0), 
F{{sq{T, fis) + 0),T, fis) < sq{T,ixb) ,V/iB > /i^- 



(12) 
(13) 



The condition (12) provides that the simple pole singularity s* = sh{T,ij,b = 0) is the 
rightmost one at vanishing /x^ = and given T, whereas the condition (13) ensures that 
s* — sq{T,ij,b) is the rightmost singularity of the isobaric partition for all values of the 
baryonic chemical potential above some positive constant /ia- This can be seen in Fig. 1 
for /ib being a variable. Since F{s,T, fiB), where it exists, is a continuous function of its 
parameters, one concludes that, if the conditions (12) and (13), are fulfilled, then at some 
chemical potential /U^(T) the both singularities should be equal. Thus, one arrives at the 
Gibbs criterion (4), but for two variables 



SH{T,fx%{T))^SQ{T,fi%{T)). 



(14) 



It is easy to see that the inequalities (12) and (13) are the sufficient conditions of a PT 
existence for more complicated functional dependencies of Fh{s,T, /ib) and Fq{s,T, /ib) 
than the ones used here. 

For the choice (9), (10) and (11) of Fh{s, T, ijlb) and -Fq(s, T, iib) functions the PT exists 
at T < Tcep, because the sufficient conditions (12) and (13) can be easily fulffiled by a proper 
choice of the bag constant B and the function u{T, hb) > for the interval T < Tup with 
the constant Tup > Tcep- Clearly, this is the 1^* order PT, since the surface tension is finite 
and it provides the convergence of the integrals (6) and (7) in the expression (5), where the 
usual T-derivatives should be now understood as the partial ones ioi hb = const. 



Assuming that the conditions (12) and (13) are fulfilled by the correct choice of the model 
parameters B and u{T, jib) > 0, one can see now that at T = Tc^p there exists a PT as well, 
but its order is defined by the value of r. As was discussed in the preceding section for 

I < r < 2 there exists the 2"*^ order PT. For 1 < r < | there exist the PT of higher order, 
defined by the conditions formulated in Eq. (8). This is a new possibility, which, to my best 
knowledge, does not contradict to any general physical principle (see the left panel in Fig. 
2). 

The case r > 2 can be ruled out because there must exist the first order PT for T > Tcep, 
whereas for T < Tcep there exists the cross-over. Thus, the critical endpoint mT — hb plane 
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Figure 2: Left panel. A schematic picture of the deconfinement phase transition diagram in 

the plane of baryonic density pb and T for the 2"'' order PT at the critical endpoint (CEP), 

i.e. for I < r < 2. For the 3'"'^ (or higher) order PT the boundary of the mixed and hadronic 

phases (dashed curve) should have the same slope as the boundary of the mixed phase and 

QGP (solid curve) at the CEP. 

Right panel. Same as in the left panel, but for r > 2. The critical endpoint in the 
fiB — T plane generates the critical end line (CELine) in the ps — T plane shown by the thick 
horizontal line. This occurs because of the discontinuity of the partial derivatives of sh and 
sq functions with respect to pb and T. 



will correspond to the critical interval in the temperature-baryonic density plane. Since such 
a structure of the phase diagram in the variables temperature-density has, to my knowledge, 
never been observed, 1 conclude that the case r > 2 is unrealistic (see the right panel in Fig. 
2). Note that a similar phase diagram exists in the EDM with the only difference that the 
boundary of the mixed and liquid phases (the latter in the QGBST model corresponds to 
QGP) is moved to infinite particle density. 



4 Surface Tension Induced Phase Transition 

Using our results for the case (III) of the preceding section, we conclude that above Tcep 
there is a cross-over, i.e. the QGP and hadrons coexist together up to the infinite values of T 
and/or pB- Now, however, it is necessary to answer the question: How can the two different 
sets of singularities that exist on two sides of the line T = Tcep provide the continuity of the 
solution of Eq. (3)? 

It is easy to answer this question for pb < p%{Tcep) because in this case all partial T 
derivatives of sh{T, Pb), which is the rightmost singularity, exist and are finite at any point of 
the line T = Tcep- This can be seen from the fact that for the considered region of parameters 
sh{T,pb) is the rightmost singularity and, consequently, sh{T,pb) > sq{T,pb). The latter 
inequality provides the existence and finiteness of the volume integral in Fq{s,T, pb)- In 
combination with the power T dependence of the reduced surface tension coefficient cr(T) 
the same inequality provides the existence and finiteness of all its partial T derivatives of 
Fq{s,T, Pb) regardless to the sign of ^(T). Thus, using the Taylor expansion in powers 
of (T — Tcep) at any point of the interval T = Tcep and pb < f^si'^cep), one can calculate 
sh{T,pb) for the values of T > Tcep which are inside the convergency radius of the Taylor 
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expansion. 

The other situation is for jis > A^s(^cep) and T > Tc^p, namely in this case above the 
deconfinement PT there must exist a weaker PT induced by the disappearance of the reduced 
surface tension coefficient. To demonstrate this we have solve Eq. (3) in the limit, when T 
approaches the curve T = T^ep from above, i.e. for T — > T^ep + 0, and study the behavior of 
T derivatives of the solution of Eq. (3) s* for fixed values of ^b- For this purpose we have to 
evaluate the integrals /Ci-(A, 7^) introduced in Eq. (7). Here the notations A = s*—sq{T, hb) 
and 7^ = — (t(T) > are introduced for convenience. 

To avoid the unpleasant behavior for r < 2 it is convenient to transform (7) further on 
by integrating by parts: 

ICr{A, f ) ^ g^iVo) - (^/C.-i(A, f ) + ^)ICr-.iA, 7^) , (15) 

where the regular function griVo) is defined as 

gr{Vo) = exp [-AVo + 7^0^ . (16) 

For T — a > 1 one can change the variable of integration v —>■ z/ A and rewrite ICr-aiA, 7^) 

I T"^ 1^ 



as 

v2 



/C._„(A, 7^) = A^-^-'l dz ^ ^ = A--«-i /C._„ (1, ^^A-^) . (17) 

VoA 

This result shows that in the limit 7^0, when the rightmost singularity must approach 
sq{T,iib) from above, i.e. A — >• 0"*", the function (17) behaves as /C^_a(A,7^) ~ ^-r-a-i _|_ 
0(A^""). This is so because for 7 ^ the ratio 7^A'"'^ cannot go to infinity, otherwise 
the function K-r-i (1,7^A~'*), which enters into the right hand side of (15), would diverge 
exponentially and this makes impossible an existence of the solution of Eq. (3) for T = T^ep- 
The analysis shows that for 7^0 there exist two possibilities: either v = 7^A~'^ Const 
or 1/ = 7^A~'^ — > 0. The most straightforward way to analyze these possibilities for 7 — > 
is to assume the following behavior 

A = A7" + 0(7"+^), ^ f = 1^ [Aar-' + 0(7")] - ^ff^jj , (18) 

and find out the a value by equating the T derivative of A with the T derivative (5). 
The analysis shows [3] that for A^""^ < 77'A-^~'' one finds 

7«-2 ^ A^-'^ ^ are = 2 for r < 1 + 3^ . (19) 

Similarly, for A^""^ > 77' A^"** one obtains 7"^~^7' ~ A^~'^ and, consequently, 

(731)^ for T>1 + ^. (20) 



Summarizing these results for 7 — > 0, one can write the expression for the second deriva- 
tive of A as [3]: 



d'A 



2k+l 



3-2t 

T-1 



r > 1 + 



2fe+l 



(21) 
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Figure 3: A schematic picture of the deconfinement phase transition diagram (full curve) 
in the plane of baryonic chemical potential /i^ and T for the 2"'' order PT at the tricritical 
endpoint (CEP). The model predicts an existence of the surface induced PT of the 2"*^ or 
higher order (depending on the model parameters). This PT starts at the CEP and goes to 
higher values of T and/or /i^. Here it is shown by the dashed curve CEP-A, if the phase 
diagram is endless, or by the dashed-dot curve CEP-B, if the phase diagram ends at T = 0, 
or by the dashed-double-dot curve CEP-C, if the phase diagram ends at /i^ = 0. Below 
(above) each of A or B curves the reduced surface tension coefficient is positive (negative). 
For the curve C the surface tension coefficient is positive outside of it. 

The last result shows us that, depending on k and k values, the second derivatives of s* and 
SQ(T,fiB) can differ from each other for | < r < 2 or can be equal for 1 < r < |. In other 
words, it is found that at the line T = Tcep there exists the 2"*^ order PT for | < r < 2 and 
the higher order PT for 1 < r < |, which separates the pure QGP phase from the region of 
a cross-over, i.e. the mixed states of hadronic and QGP bags. Since it exists at the line of a 
zero surface tension, this PT will be called the surface induced PT. For instance, from (21) 
it follows that for /c = and > | there is the 2""^ order PT, whereas for = and = | 
or for k > and k, < 1 there is the 3'^ order PT, and so on. 

Since the analysis performed in the present section did not include any /i^ derivatives of 
A, it remains valid for the hb dependence of the reduced surface tension coefficient, i.e. for 
Tcep{p-B)- Only it is necessary to make a few comments on a possible location of the surface 
tension null line Tcepi^^s) ■ In principle, such a null line can be located anywhere, if its location 
does not contradict to the sufficient conditions (12) and (13) of the I''* deconfinement PT 
existence. Thus, the surface tension null line must cross the deconfinement line in the ^s—T 
plane at a single point which is the tricritical endpoint (/x^^; Tcep(/i^^)), whereas for /i^ > /^^^ 
the null line should have higher temperature for the same than the deconfinement one, 
i.e. TcepifJ^B) > Tc^Hb) (see Fig. 3). Clearly, there exist two distinct cases for the surface 
tension null line: either it is endless, or it ends at zero temperature or at other singularity, 
like the Color-Flavor-Locked phase. From the present lattice QCD data [32] it follows that 
the case C in Fig. 3 is the least possible. 

To understand the meaning of the surface induced PT it is instructive to quantify the 
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difference between phases by looking into the mean size of the bag: 



(22) 

s=s*-0 

As was shown in hadronic phase phase A > and, hence, it consists of the bags of finite 
mean volumes, whereas, by construction, the QGP phase is a single infinite bag. For the 
cross-over states A > and, therefore, they are the bags of finite mean volumes, which 
gradually increase, if the rightmost singularity approaches sq{T,ij,b), i-e. at very large 
values T and/or fiB- Such a classification is useful to distinguish QCD phases of present 
model: it shows that hadronic and cross-over states are separated from the QGP phase by 
the 1** order deconfinement FT and by the 2"°' or higher order FT, respectively. 

5 Concluding Remarks 

Here I presented an analytically solvable statistical model which simultaneously describes 
the 1** and 2"*^ order PTs with a cross-over. The approach is general and can be used for 
more complicated parameterizations of the hadronic mass- volume spectrum, if in the vicinity 
of the deconfinement FT region the discrete and continuous parts of this spectrum can be 
expressed in the form of Eqs. (9) and (10), respectively. Also the actual parameterization 
of the QGP pressure p — Tsq{T, hb) was not used so far, which means that our result can 
be extended to more complicated functions, that can contain other phase transformations 
(chiral FT, or the FT to color superconducting phase) provided that the sufficient conditions 
(12) and (13) for the deconfinement FT existence are satisfied. 

In this model the desired properties of the deconfinement phase diagram are achieved 
by accounting for the temperature dependent surface tension of the quark-gluon bags. As 
was shown, it is crucial for the cross-over existence that at T = T^ep the reduced surface 
tension coefficient vanishes and remains negative for temperatures above T^ep- Then the 
deconfinement jiB — T phase diagram has the FT at hb > f^BC^cep) for | < r < 2 , which 
degenerates into the 2"^^ order FT (or higher order FT for | > r > 1) at /^b = ii%{Tcep), and 
a cross-over for < /^b < ii%{Tcep)- These two ingredients drastically change the critical 
properties of the GBM [17] and resolve the long standing problem of a unified description 
of the 1'* and T"^ order FTs and a cross-over, which, despite all claims, was not resolved in 
Ref. [34]. In addition, it was found that at the null line of the surface tension there must 
exist the surface induced FT of the 2"^^ or higher order, which separates the pure QGF from 
the mixed states of hadrons and QGF bags, that coexist above the cross-over region (see 
Fig. 3). Thus, the QGBST model predicts that the QCD critical endpoint is the tricritical 
endpoint. It would be interesting to verify this prediction with the help of the lattice QCD 
analysis. For this one will need to study the behavior of the bulk and surface contributions 
to the free energy of the QGF bags and/or the string connecting the static quark-antiquark 
pair. 

However, the QGP bags created in the experiments have finite mass, volume, life-time 
and, hence, the strong discontinuities which are typical for the P* order FT should be 
smeared out which would make them hardly distinguishable from the cross-over. Thus, to 
seriously discuss the signals of the F** order deconfinement FT and/or the tricritical endpoint. 



^ ainF(s,T,MB) 
' d s 



(v) 
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one needs to solve the finite volume version of the QGBST model hke it was done for the 
SMM [11] and the GBM [12, 35]. This, however, is not sufficient because, in order to make 
any reliable prediction for experiments, the finite volume equation of state must be used 
in hydrodynamic equations which, unfortunately, are not suited for such a purpose. Thus, 
we are facing a necessity to return to the foundations of heavy ion phenomenology and to 
modify them according to the requirements of the experiments. 

In addition, to apply the QGBST model to the experiments it is nesseary to make it 
more realistic: it seems that for the mixture of hadrons and QGP bags above the cross-over 
line it is necessary to include the relativistic treatment of hard core repulsion [36, 37] for 
lightest hardons and to include into statistical description the medium dependent width of 
hadronic resonances and QGP bags, which, as argued in Ref. [38], may completely change 
our understanding of the cross-over mechanism. 
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